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Abstract. Given a set of n terminals, which are points in d-dimensional Euchdean space, the mmimum 
Manhattan network problem (MMN) asks for a minimum-length rectilinear network that connects each 
pair of terminals by a Manhattan path, that is, a path consisting of axis-parallel segments whose total 
length equals the pair's Manhattan distance. Even for d — 2, the problem is NP-hard, but constant- 
factor approximations are known. For d > 3, the problem is APX-hard; it is known to admit, for any 
e > 0, an 0(n'^ )-approximation. 

In the generalized minimum Manhattan network problem (GMMN), we are given a set R oi n terminal 
pairs, and the goal is to find a minimum-length rectilinear network such that each pair in R is connected 
by a Manhattan path. GMMN is a generalization of both MMN and the well-known rectilinear Steiner 
arborescence problem (RSA). So far, only special cases of GMMN have been considered. 
We present an 0(log'''''^ n)-approximation algorithm for GMMN (and, hence, MMN) in d > 2 dimen- 
sions and an 0(log n)-approximation algorithm for 2D. We show that an existing 0(log n)-approximation 
algorithm for RSA in 2D generalizes easily to d > 2 dimensions. 



1 Introduction 

Given a set of terminals, which are points in M , the minimum Manhattan network problem (MMN) 
asks for a minimum-length rectilinear network that connects every pair of terminals by a Manhattan 
path (M-path, for short), that is, a path consisting of axis-parallel segments whose total length equals 
the pair's Manhattan distance. 

In the generalized minimum Manhattan network problem 
(GMMN), we are given a set Roin unordered terminal pairs, 
and the goal is to find a minimum-length rectilinear network 
such that every pair in R is M-connected, that is, connected 
by an M-path. GMMN is a generalization of MMN since R 
may contain all possible pairs of terminals. Figure [T] depicts 
examples of both network types. 

We remark that in this paper we define n to be the num- 
ber of terminal pairs of a GMMN instance, whereas previous 
works on MMN defined n to be the number of terminals. 

Two-dimensional MMN (2D-MMN) naturally arises in 
VLSI circuit layout [GLNOlj . where a set of terminals (such as gates or transistors) needs to be 
interconnected by rectilinear paths (wires) . Minimizing the cost of the network (which means min- 
imizing the total wire length) is desirable in terms of energy consumption and signal interference. 
The additional requirement that the terminal pairs are connected by shortest rectilinear paths aims 
at decreasing the interconnection delay (see Cong et al. |CLZ93j for a discussion in the context of 
rectilinear Steiner arborescences, which have the same additional requirement; see definition below). 







f 


7 
d 


□« 




c 








,b 











y 

o" 


/ 


Co 

c 





(a) an MMN for 
{a,b,c,d,e,f} 



(b) a GMMN for 
{(a,fe),(c,d),(e,/)} 



Fig. 1: MMN versus GMMN in 2D. 



Manhattan networks also arise in the area of geometric spanner networks. Specifically, a minimum 
Manhattan network can be thought of as the cheapest spanner under the Li-norm for given a set of 
points (allowing Steiner points). Spanners, in turn, have numerous applications in network design, 
distributed algorithms, and approximation algorithms. 

MMN requires a Manhattan path between every terminal pair. This assumption is, however, 
not always reasonable. Specifically in VLSI design a wire connection is necessary only for a, often 
comparatively small, subset of terminal pairs, which may allow for substantially cheaper circuit 
layouts. In this scenario, GMMN appears to be a more realistic model than MMN. 

Previous Work. MMN was introduced by Gudmundsson et al. |GLN01| who gave 4- and 8- 
approximation algorithms for 2D-MMN running in 0{n^) and 0(n log n) time, respectively. The 
currently best known approximation algorithms for 2D-MMN have ratio 2; they were obtained 
independently by Chepoi et al. |CNV08| using an LP-based method, by Nouioua |Nou05] using a 
primal-dual scheme, and by Guo et al. |GSZ11) using a greedy approach. The complexity of 2D- 
MMN was settled only recently by Chin et al. |CGS11| : they proved the problem NP-hard. It is 
not known whether 2D-MMN is APX-hard. 

Less is known about MMN in dimensions greater than 2. Muhoz et al. JMSU09J proved that 
3D-MMN is NP-hard to approximate within a factor of 1.00002. They also gave a constant-factor 
approximation algorithm for a, rather restricted, special case of 3D-MMN. More recently. Das et 
al. described the first approximation algorithm for MMN in arbitrary, fixed dimension with a ratio 
of Oin") for any e > |DGK+llj . 

GMMN was defined by Chepoi et al. jCNV OB] who asked whether 2D-GMMN admits an 0(1)- 
approximation. Apart from the formulation of this open problem, only special cases of GMMN 
(such as MMN) have been considered so far. 

One such special case (other than MMN) that has received significant attention in the past is 
the rectilinear Steiner arborescence problem (RSA). Here, we are given n terminals lying in the 
first quadrant and the goal is to find a minimum-length rectilinear network that M-connects every 
terminal to the origin o. Hence, RSA is the special case of GMMN where o is considered a (new) 
terminal and the set of terminal pairs contains, for each terminal t ^ o, only the pair (o, t). RSA 
was introduced by Nastansky et al. |NSS74j and has mainly been studied in 2D. 2D-RSA is NP- 
hard jSSOO] . Rao et al. |RSHS92] gave a 2-approximation algorithm for 2D-RSA. They also provided 
a conceptually simpler 0(logn)-approximation algorithm based on rectilinear Steiner trees. That 
algorithm generalizes quite easily to dimensions d > 2 (as we show in the appendix). Lu et al. |LR00j 
and, independently, Zachariasen [ZacOOj described polynomial-time approximation schemes (PTAS) 
for RSA, both based on Arora's technique [ Aro03] . Zachariasen pointed out that his PTAS can be 
generalized to the all-quadrant version of RSA but that it seems difficult to extend the approach 
to higher dimensions. 

Our Contribution. Our main result is an 0(log ^ n)-approximation algorithm for GMMN (and, 
hence, MMN) in d dimensions. For the sake of simplicity, we first present our approach in 2D (see 
Section [2]) and then show how it can be generalized to higher dimensions (see Section^. We also 
provide an improved and technically more involved 0(logn)-approximation for the special case of 
2D-GMMN, but this approach does not seem to generalize to higher dimensions; see Section |4] To 
the best of our knowledge, we present the first approximation algorithms for GMMN. Our result 
for 2D is not quite the constant-factor approximation that Chepoi et al. were asking for, but it 
is a considerable step into that direction. Note that the poly-logarithmic ratio of our algorithm 



for GMMN in d > 3 dimensions constitutes an exponential improvement upon the previously only 
known approximation algorithm, which solves the special case MMN, with a ratio of 0{rf) for any 



e >0 DGK+11 



Our algorithm for GMMN is based on divide and conquer. We identify each terminal pair with 
its d-dimensional bounding box. Consequently, we consider R a set of d-dimensional boxes. We use 
(d — l)-dimensional hyperplanes to partition R recursively into sub- instances. The base case of our 
partition scheme consists of GMMN instances where all boxes contain a common point. We solve 
the resulting special case of GMMN by reducing it to RSA. We have postponed the running-time 
analysis to Appendix [B| 

2 Polylogarithmic Approximation for Two Dimensions 

In this section, we present an 0(log n)-approximation algorithm for 2D-GMMN. The algorithm 
consists of a main algorithm that recursively subdivides the input instance into instances of so- 



called x-separated GMMN; see Section 2.1 We prove that the instances of x-separated GMMN 



can be solved independently by paying a factor of O(logn) in the approximation ratio. Then we 



show how to approximate x-separated GMMN within ratio O(logn); see Section 2.2 This yields 
an overall ratio of 0(log n). 

2.1 Main Algorithm 

Our approximation algorithm is based on divide and conquer. Let R be the set of terminal pairs 
that are to be M-connected. We identify each terminal pair with its bounding box, that is, the 
smallest axis-aligned rectangle that contains both terminals. As a consequence of this, we consider 
R a set of rectangles. Let m^ be the median in the multiset of the x-coordinates of terminals. We 
identify rrix with the vertical line at x = nix- 

Our algorithm divides R into three subsets i?ieft , -Rmid , and -Rright • The set i?ieft consists of all 
rectangles that lie completely to the left of the vertical line nix- Similarly, the set i?right consists 
of all rectangle that lie com,pletely to the right of rux- The set -Rmid consists of all rectangles that 
intersect rrix. 

We consider the sets -Rieft, Rmid-, and -Rright as separate instances of GMMN and apply our 
algorithm recursively to R\e^ and to -Rright • The union of the two resulting networks is a rectilinear 
network that M-connects all terminal pairs in i?icft U -Rright • 

It remains to M-connect the pairs in -Rmid- We call an GMMN instance (such as i?mid) x- 
separated if there is a vertical line (in our case rUx) that intersects every rectangle. We exploit 
this property to design a simple 0(logn)-approximation algorithm for x-separated GMMN; see 
Section 2.2 Later, in SectionHl we improve upon this and describe an 0(l)-approximation algorithm 



for x-separated GMMN. 

To analyze the performance of our main algorithm, let p{n) denote the algorithm's worst-case 
approximation ratio for instances with n terminal pairs. Now assume that our input instance R is 
a worst case. More precisely, the cost of the solution of our algorithm equals p{n) ■ OPT, where 
OPT denotes the cost of an optimum solution A^opt to R. Let Aioft and N^ight be the parts of Aopt 
to the left and to the right of nix, respectively. 

Due to the choice of rUx, at most n terminals lie to the left of rUx- Therefore, -Rieft contains at 
most n/2 terminal pairs. Since Aieft is a feasible solution to i?icft, we conclude that the cost of the 
solution to i?icft computed by our algorithm is bounded by p{n/2) ■ ||Aicft||, where || • || measures 



the length of a network. Analogously, the cost of the solution computed for Rnght is bounded by 
p{n/2) ■ llA^^rightll- Now we assume that we can approximate x-separated instances with a ratio of 
Px{n). Since A'^Qpt is also a feasible solution to the x-separated instance -Rmid; we can compute a 
solution of cost px{n) ■ OPT for .Rmid- 

Therefore, we can bound the total cost of our algorithm's solution A^ to i? by 

p{n) • OPT = ||iV|| < p(n/2) • (||Afieft|| + ||iVright||) + Px{n) ■ OPT . 

Note that this inequality does not necessarily hold if R is not a worst case since then p(n) ■ 
OPT > ||A^||. The networks Aieft and Alright are separated by nix, hence they are edge disjoint and 
||A^ieft|| + llA^rightll < OPT. This yields the recurrence p{n) < p{n/2) + px{n), which resolves to 
p{n) = logn ■ Px{n). Let's summarize this discussion. 

Lemma 1. If x- separated 2D-GMMN admits a px{n)- approximation, 2D-GMMN admits a {px{n) ■ 
log n)- approximation. 

Combining this lemma with our O(logn) approximation algorithm for x-separated instances 
described below, we obtain the following intermediate result. 

Theorem 1. 2D-GMMN admits an 0(log n)- approximation. 

2.2 Approximating a;-Separated Instances 

In this section, we describe a simple algorithm for approximating x-separated 2D-GMMN instances 
with a ratio of O(logn). Let R be our input. Since R is x-separated, all rectangles in R intersect a 
common vertical line. W.l.o.g., this is the y-axis. 

The algorithm works as follows. Analogously to the main algorithm presented in Section 2.1 , we 
recursively subdivide the x-separated input instance, but this time according to the y- coordinate. 
As a result of this, the input instance R is decomposed into y-separated sub-instances. Moreover, 
since each of these sub-instances is (as a subset of R) already x-separated, we call these instances 



xy-separated. In Section [ZSj we give a specialized algorithm for xy-separated instances. 

Let Px{n) be the ratio of our algorithm for approximating x-separated GMMN instances and 
let Pxy{n) be the ratio of our algorithm for approximating xy-separated GMMN instances. In 



Section 2.3, we show that Pxy{n) = 0(1)- Then Lemma IT] (by exchanging x- and y-coordinates) 



implies that Px{n) = logn ■ Pxy{n) = O(logn). 

Lemma 2. x-separated 2D-GMMN admits an 0{\ogn)-approximation. 

2.3 Approximating a;y-Separated Instances 

It remains to show that xy-separated GMMN instances can be approximated within a constant 
ratio. Let R be such an instance. We assume, w.l.o.g., that it is the x- and the y-axis that intersect 
all rectangles in i?, that is, all rectangles contain the origin. Let Aopt be an optimum solution to R. 
Let A^ be the union of Aopt with the projections of A^opt to the x-axis and to the y-axis. The total 
length of A^ is ||A^|| < 2 • OPT = O(OPT) since every line segment of Aopt is projected either to 
the X-axis or to the y-axis but not to both. The crucial fact about A^ is that this network contains, 
for every terminal t in R, an M-path from t to the origin a. In other words, A^ is a feasible solution 
to the RSA instance of M-connecting every terminal in R to o. 



To see this, consider an arbitrary terminal pair {t,t') G R. Let 77 be yk 
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an M-path connecting t and t in A^opt; see Fig. YA Note that, since the 

bounding box of (t, t') contains o, 77 intersects both x- and y-axis. To 

obtain an M-path from t to o, we follow 77 from t to t' until 77 crosses 

one of the axes. From that point on, we follow the projection of 77 on 

this axis. We reach o when 77 crosses the other axis; see the dotted path i* 

in Fig. [2] Analogously, we obtain an M-path from t' to 0. Yis,. 2: Network A'' connects 

Let T be the set of terminals in R. We have shown above that there is all terminals to origin o. 
a feasible solution N of cost O(OPT) to the RSA instance with terminal 

set T. There is a PTAS for RSA in two dimensions [ZacOOl ILROO] . Using this PTAS, we can 
efficiently compute a feasible RSA solution A^' for T of cost 0(1) • ||Ar|| = O(OPT). Moreover, A^' 
is also a feasible solution to the GMMN instance 7?. To see this, note that A^' contains, for every 
terminal pair {t,t') G 7?, an M-path vr from t to o and an M-path vr' from to t' . Concatenating tt 
and vr' yields an M-path from t to t' as the bounding box of (t, t') contains o. Thus we obtain the 
following result. 

Lemma 3. xy-separated 2D-GMMN admits a constant-factor approximation. 

3 Generalization to Higher Dimensions 

In this section, we describe an 0(log ^ n)-approximation algorithm for GMMN in d dimensions, 
which is a generalization of the algorithm for two dimensions presented in Section [2| Let us view 



this algorithm from the following perspective. In Section 2.1, we reduced GMMN to solving x- 



separated sub-instances at the expense of a (log n)-factor in the approximation ratio (see Lemma [I|. 



Applying the same lemma to the y-coordinates in Section 2.2, we further reduced the problem to 
solving j;y-separated sub-instances, that is, to instances that were separated with respect to both 
dimensions. This caused the second (log n)-factor in our approximation ratio. Finally, we were able 
to approximate these completely separated sub-instances within constant ratio by solving a related 



RSA problem (see Section 2.3) 



These ideas generalize to higher dimensions. An instance 7? of d-dimensional GMMN is called 
j -separated for some j < d if there exist values si, . . . ,Sj such that, for each terminal pair (t, t') £ R 
and for each dimension i < j, we have that Sj separates the i-th coordinates Xi{t) of t and Xj(t') 
of t' (meaning that either Xi{t) < Sj < Xj(t') or Xi{t') < Sj < Xi(t)). Under this terminology, an 
arbitrary instance of d-dimensional GMMN is always 0-separated. 

We first show that if we can approximate j-separated GMMN with ratio Pj{n) then we can 



approximate {j — l)-separated GMMN with ratio Pj{n) ■ logn; see Section 3.1 Then we show 
that d-separated GMMN can be approximated within a factor pd{n) = O(logn); see Section 3.2 
Combining these two facts and applying them inductively to arbitrary (that is, 0-separated) GMMN 
instances yields the following central result of our paper. 

Theorem 2. GMMN in d dimensions admits an 0(log ^ n) -approximation. 

As a byproduct of this algorithm, we obtain an 0(log "*" n)-approximation algorithm for MMN 
where n denotes the number of terminals. This holds since any MMN instance with n terminals 
can be considered an instance of GMMN with O(n^) terminal pairs. 

Corollary 1. MMN in d dimensions admits an 0(log ^ n)- approximation, where n denotes the 
number of terminals. 



3.1 Separation 

In this section, we show that if we can approximate j-separated GMMN instances with ratio Pj{n), 
we can approximate {j — l)-separated instances with ratio logn • Pj{n). The separation algorithm 
and its analysis work analogously to the main algorithm for 2D where we reduced (approximating) 



2D-GMMN to (approximating) x-separated 2D-GMMN; see Section [23 

Let i? be a set of {j — l)-separated terminal pairs. Let rrix be the median in the multiset of the 
j-th. coordinates of terminals. We divide R into three subsets i?icft, ^mid) and i?right- The set i?ieft 
consists of all terminal pairs {t,t') such that Xj{t),Xj{t') < nix and -Rright contains all terminal 
pairs (i, t') with Xj{t),Xj{t') > rrix. The set -Rmid contains the remaining terminal pairs, all of which 
are separated by the hyperplane Xj = rux- We apply our algorithm recursively to Rmi and i?right- 
The union of the resulting networks is a rectilinear network that M-connects all terminal pairs 

Rleft U -Rright- 

In order to M-connect the pairs in -Rmid, we apply an approximation algorithm for j-separated 
GMMN of ratio Pj{n). Note that the instance -Rmid is in fact j-separated by construction. The 
analysis of the resulting algorithm for {j — l)-separated GMMN is analogous to the 2D-case (see 



Section 2.1) and is therefore omitted. 



Theorem 3. Let I < j < d. If j -separated GMMN admits a pj{ri)- approximation, then {j — 1)- 
separated GMMN admits a {pj{n) ■ logn) -approximation. 

3.2 Approximating d-Separated Instances 

In this section, we show that we can approximate instances of d-separated GMMN within a ratio of 
O(logn) by reducing the problem to RSA. Let -R be a d-separated instance and let T be the set of 
all terminals in R. As R is d-separated, all bounding boxes defined by terminal pairs in R contain 
a common point, which is, w.l.o.g., the origin. 



As in the two-dimensional case (see Section 2.3), we M-connect all terminals to the origin by 



solving an RSA instance with terminal set T. This yields a feasible GMMN solution to R since for 
each pair {t,t') G R there is an M-path from t to the origin and an M-path from the origin to t' . 
The union of these paths is an M-path from t to t' since the origin is contained in the bounding 
boxof (i,t'). 

Rao et al. |RSHS92] presented an 0(log |T|)-approximation algorithm for 2D-RSA, which gener- 
alizes, in a straight-forward manner, to d-dimensional RSA; see Appendix [A] for details. Hence, we 
can use the algorithm of Rao et al. to efficiently compute a feasible GMMN solution. The following 
lemma shows that this solution is in fact an 0(logn)-approximation. The proof is similar to the 
proof of Lemma 7 in the paper of Das et al. |DGK"'"ll] . 

Lemma 4. d-separated GMMN admits an O {log n)- approximation for any fixed number d of di- 
mensions. 

Proof. Below we show that there is a solution of cost O(OPT) to the RSA instance connecting T 
to the origin. Observing that \T\ < 2n and using our extension of the 0(log |T|)-approximation 
algorithm of Rao et al. (see Appendix [A|), we can efficiently compute a feasible GMMN solution of 
cost O(OPT-logn). 

We now show that there is an RSA solution of cost O(OPT). Let A'opt be an optimal GMMN 
solution to R and let A^ be the projection of A'opt onto all subspaces that are spanned by some 

6 



subset of the coordinate axes. Since there are 2*^ such subspaces, which is a constant for fixed d, 
the cost of A^is O(OPT). 

It remains to show that N M-connects all terminals to the origin, that is, A^ is a feasible solution 
to the RSA instance. First, note that A^opt ^ A^ since we project on the d-dimensional space, too. 
Now consider an arbitrary terminal pair (t, t') in R and an M-path vr in A^opt that M-connects t and 
t' . Starting at t, we follow vr until we reach the first point pi where one of the coordinates becomes 
zero. W.l.o.g., xi{pi) = 0. Clearly vr contains such a point as the bounding box of (t,t') contains 
the origin. We observe that pi lies in the subspace spanned by the d — 1 coordinate axes X2, ■ ■ ■ , Xd- 
From pi on we follow the projection of vr onto this subspace until we reach the first point p2 where 
another coordinate becomes zero; w.l.o.g., X2{p2) = 0. Hence, p2 has at least two coordinates that 
are zero, that is, p2 lies in a subspace spanned by only d — 2 coordinate axes. Iteratively, we continue 
following the projections of vr onto subspaces with a decreasing number of dimensions until every 
coordinate is zero, that is, we have reached the origin. An analogous argument shows that A^ also 
contains an M-path from t' to the origin. D 

4 Improved Algorithm for Two Dimensions 

In this section, we show that 2D-GMMN admits an 0(logn)-approximation, which improves upon 
the 0(log n)-result of Sectional To this end, we develop a (6 + e)-approximation algorithm for 
x-separated 2D-GMMN, for any e > 0. While the algorithm is simple, its analysis turns out to be 
quite intricate. In Appendix [Cj we show tightness. Using Lemma [T| our new subroutine for the 
j;-separated case yields the following. 

Theorem 4. 2D-GMMN admits a ((6 + e) ■ log n)- approximation. 

Let R be the set of terminal pairs of an x-separated instance of 2D-GMMN. We assume, w.l.o.g., 
that each terminal pair {l,r) £ R is separated by the y-axis, that is, x{l) < < x(r). Let A'opt be 
an optimum solution to R. Let OPTver and OPThor be the total costs of the vertical and horizontal 
segments in A'opt, respectively. Hence, OPT = OPTver + OPThor- We first compute a set S of 
horizontal line segments of total cost O(OPThor) such that each rectangle in R is stabbed by some 



line segment in S; see Sections 4.1 and|4.2[ Then we M-connect the terminals to the y-axis so that 



the resulting network, along with the affected part of the y-axis and the stabbing S, forms a feasible 



solution to R of cost O(OPT); see Section 4.3 



4.1 Stabbing the Right Part 

We say that a horizontal line segment h stabs an axis-aligned rectangle r if /i intersects the boundary 
of r twice. A set of horizontal line segments is a stabbing of a set of axis-aligned rectangles if each 
rectangle is stabbed by some line segment. For any geometric object, let its right part be its 
intersection with the closed half plane to the right of the y-axis. For a set of objects, let its right 
part be the set of the right parts of the objects. Let i?"*" be the right part of R, let A^"*" be the right 
part of A'opt, and let A^j^^. be the set of horizontal line segments in A^^. In this section, we show 
how to construct a stabbing of R'^ of cost at most 2 • IJA^j^^H. 

For x' > 0, let i^' be the vertical line at x = x' . Our algorithm performs a left-to-right sweep 
starting with Iq. Note that, for every x > 0, the intersection of R'^ with Ix forms a set Ix of 
intervals. The intersection of A^^^^. (|J N^^^ to be precise) with ix is, at any time, a set of points 



that constitutes a piercing for I^, that is, every interval in I^ contains a point in ix Pi N^^^. Note 
that||iV+J|=/|4niV+Jdx. 

We imagine that we continuously move ix from x = to the right. At any time, we maintain 
an inclusion-wise minimal piercing Px oi Xx- With increasing x, we only remove points from Px\ 
we never add points. This ensures that the traces of the points in Px form horizontal line segments 
that all touch the y-axis. These line segments form our stabbing of R^ . 

The algorithm proceeds as follows. It starts at x := with an arbitrary minimal piercing Pq 
of Xq. Note that we can even compute an optimum piercing Pq. We must adapt Px whenever Xx 
changes. With increasing x, Xx decreases inclusion-wise since all rectangles in R^ touch the y-axis. 
So it suffices to adapt the piercing Px at event points; x is an event point if and only if x is the 
x-coordinate of a right edge of a rectangle in i?"*". 

Let x' and x" be consecutive event points. Let x be such that x' < x < x" . Note that Px' is 
a piercing for Xx since Xx C Xx'- The piercing Px' is, however, not necessarily minimal w.r.t. Xx- 
When the sweep line passes x', we therefore have to drop some of the points in P^./ in order to 
obtain a new minimal piercing. This can be done by iteratively removing points from Px' such that 
the resulting set still pierces Xx- We stop at the last event point (afterwards, Xx = 0) and output 
the traces of the piercing. 

It is clear that the algorithm produces a stabbing of R^] see the thick solid line segments in 



Fig. 3a The following lemma is crucial to prove the overall cost of the stabbing. 



Lemma 5. For any x >0, it holds that \Px\ < 2 • |£^ n -/Vj^j,|. 

Proof. Since Px is a minimal piercing, there exists, for every p £ Px, a witness interval Ip G Xx 
that is pierced by p but not by Px \ {p} ■ Otherwise we could remove p from Px , contradicting the 
minimality of Px- 

Now we show that an arbitrary point q on ix is contained in the witness intervals of at most 
two points in Px- Assume, for the sake of contradiction, that q is contained in the witness intervals 
of points p,p',p" G Px with strictly increasing y-coordinates. Suppose that q lies above p' . But 
then the witness interval Ip of p, which contains p and q, must also contain p' , contradicting the 
definition of Ip. The case q below p' is symmetric. 

Recall that ix H N^^^ is a piercing of Xx and, hence, of the \Px\ many witness intervals. Since 
every point in ix Pi -/Vj^j, pierces at most two witness intervals, the lemma follows. D 

Observe that the cost of the stabbing is J \Px\dx. By the above lemma, the cost of the stabbing 
can be bounded by f \Px\dx < f 2 ■ \ixf^ -^horl*^-^ ~ ^ " H^^horll' '^hich proves the following lemma. 

Lemma 6. Given a set R of rectangles intersecting the y-axis, we can compute a set of horizontal 
line segments of cost at most 2 • OPThor that stabs R~^ . 



4.2 Stabbing Both Parts 

We now detail how we construct a stabbing of R. To this end we apply Lemma |6] to compute a 
stabbing S~ of cost at most 2 • ||A^ho;,|| for the left part R~ of R and a stabbing S"^ of cost at most 
2 • llA'^j^j.ll for the right part R~^. Note that S~ U S'^ is not necessarily a stabbing for R since there 
can be rectangles that are not completely stabbed by one segment. To overcome this difficulty, we 



mirror S and 5""*" to the respective other side of the y-axis; see Fig. 3a The total cost of the 
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Fig. 3: The improved algorithm for a;-separated 2D-GMMN. 



resulting set S of horizontal line segments is at most 4(||A'j^^|| + IJA'^^j.ll) = 4 • OPThor- The set S 
stabs R since, for every rectangle r £ R, the larger among its two (left and right) parts is stabbed 
by some segment s and the smaller part is stabbed by the mirror image s' of s. Hence, r is stabbed 
by the line segment sL) s' . Let us summarize. 

Lemma 7. Given a set R of rectangles intersecting the y-axis, we can compute a set of horizontal 
line segments of cost at most 4 • OPThoj. that stabs R. 



4.3 Connecting Terminals and Stabbing 

We assume that the union of the rectangles in R is connected. Otherwise we apply our algorithm 
separately to each subset of R that induces a connected component of \J R. Let I be the line 
segment that is the intersection of the y-axis with |J .R. Let top(/) and bot(/) be the top and 
bottom endpoints of /, respectively. Let L be the set containing every terminal t with (t, t') £ R 
and y{t) < y{t'). Symmetrically, let H be the set containing every terminal t with {t,t') G R and 
y{t) > y{t'). Note that L and H are not necessarily disjoint. 

Using a PTAS for RSA [LROOl IZacOO] , we compute a near-optimal RSA network A^p connecting 
the terminals in L to top (I) and a near-optimal RSA network j4down connecting the terminals in H 
to bot(/). Then we return the network N = ^up U ^down U S, where S is the stabbing computed 



by the algorithm in Section 4.2 



We now show that this network is a feasible solution to R. Let (/, h) E R. W.l.o.g., / G L and 
h G H. Hence, A^p contains a path tti from / to top(/), see Fig. 3b, This path starts inside the 



rectangle {l,h) and leaves it through its top edge. Before leaving {l,h), the path intersects a line 
segment s in 5 that stabs (/, h). This line segment is also intersected by the path n^ in ^down that 
connects h to bot(I). Hence, walking along tti, s, and iTh brings us in a monotone fashion from / 
to h. 

Now, let us analyze the cost of A^. Clearly, the projection of Aopt onto the y-axis yields the line 
segment I. Hence, |/| < OPTver- Observe that AoptU{/} constitutes a solution to the RSA instance 
(L, top(/)) connecting all terminals in L to top(I) and to the RSA instance {H,hot{I)) connecting 
all terminals in H to bot(/). This holds since, for each terminal pair, its M-path vr in Aopt crosses 
the y-axis in /; see Fig. 3c Since A^p and ^down are near-optimal solutions to these RSA instances, 
we obtain, for any e > 0, that \\Anp\\ < (1 + e) • ||Aopt U /|| < (1 + e) • (OPT + OPTver) and 
analogously pdownll < (1 + e) • (OPT + OPTver)- 



By Lemma ^ we have that ||5|| < 4 • OPThor- Assuming e < 1, this yields 

1 1 iV 1 1 = II ^up 11 + 11 ^down II ~r 1 1 '-' 1 1 

< (2 + 2e) ■ (OPT + OPTver) + 4 • OPThor 

< (2 + 2e) ■ OPT + 4 • (OPT^cr + OPThor) 
= (6 + e') ■ OPT 

for e' = e/2, which we can make arbitrarily small by making e arbitrarily small. We summarize our 
result as follows. 

Lemma 8. x-separated 2D-GMMN admits, for any e > 0, a (6 + e) -approximation. 

5 Conclusion and Open Problems 

We have presented an 0(log ^ n)-approximation algorithm for d-dimensional GMMN, which im- 
plies the same ratio for MMN. Prior to our work, no approximation algorithm for GMMN was 
known. For d > 3, our result is a significant improvement over the ratio of O(n^) of the only 
approximation algorithm for d-dimensional MMN known so far. 

In 2D, there is still quite a large gap between the currently best approximation ratios for MMN 
and GMMN. Whereas we have presented an 0(logn)-approximation algorithm for 2D-GMMN, 
2D-MMN admits 2- approximations [(^NVOSl IGSZlll INoun5j — but is 2D-GMMN really harder to 
approximate than 2D-MMN? Indeed, given that GMMN is more general than MMN, it may be 
possible to derive stronger non-approximability results for GMMN. So far, the only such result is 
that 3D-MMN cannot be approximated beyond a factor of 1.00002 |MSU09j . 

Concerning the positive side, for d > 3, a constant-factor approximation for d-dimensional RSA 
would shave off a factor of O(logn) from the current ratio for d-dimensional GMMN. This may be 
in reach given that 2D-RS A admits even a PTAS |LR001 IZacOO] . Alternatively, a constant-factor 
approximation for [d — A;)-separated GMMN for some k < d would shave off a factor of 0(log n) 
from the current ratio for d-dimensional GMMN. 
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Appendix 

A Solving RSA in Higher Dimensions 

In this section, we show that we can approximate d-dimensional RSA with a ratio of O(logn) even 
in the ah-orthant case where every orthant may contain terminals. In this section, n denotes the 
number of terminals. We generahze the algorithm of Rao et al. |RSHS92 ] who give an O(logn)- 
approximation algorithm for the one-quadrant version of 2D-RSA. 

It is not hard to verify that the 0(logn)-approximation algorithm of Rao et al. carries over 
to higher dimensions in a straightforward manner if all terminals lie in the same orthant. We 
can therefore obtain a feasible solution to the all-orthant version by applying the approximation 
algorithm to each orthant separately. This worsens the approximation ratio by a factor no larger 
than 2°^ since there are 2'^ orthants. Hence, we can quite easily give an 0(logn)-approximation 
algorithm for the all-orthant version since 2"^ is a constant for fixed dimension d. 

In what follows, we present a tailor-made approximation algorithm for the all-orthant version 
of d-dimensional RSA that avoids the additional factor of 2"^. Our algorithm is an adaption of the 
algorithm of Rao et al., and our presentation closely follows their lines, too. 

Consider an instance of RSA given by a set T of terminals in M'^ (without restriction of the 
orthant). Let o denote the origin. The algorithm relies on the following lemma, which we prove 
below. 

Lemma 9. Given a rectilinear Steiner tree B for terminal set TU{o}, we can find an RSA network 
A for T of length at most [log2 ri\ ■ \\B\\. 

Every RSA network is also a rectilinear Steiner tree. Since the rectilinear Steiner tree problem 
(RST) admits a PTAS for any fixed dimension d |Aro98j . we can generate a (1 + e)-approximate 
RST network B that connects T and the origin. By means of Lemma ^ we get a {1 + e) [log2 n\- 
approximation for the RSA instance T. 

Theorem 5. The all-orthant version of d- dimensional RSA admits a (l+e)- [log2 n] approximation 
for any e > 0. 

Our proof of Lemma [9] relies on the following technical lemma, which constitutes the main 
modification that we make to the algorithm of Rao et al. See Fig. |4] for an illustration. 

Lemma 10. Let t, t' be two terminals. Then we can compute in constant time a point min(f , t') 
and an M-path Tr{t, t') from t to t' containing min(t, t') with the following property. The union of 
7r(t, t') with an M-path from min(t, t') to o M-connects t and t' to o. 

Proof. We start with the following simple observation. If s and s' are points and p is a point in the 
bounding box -B(s, s') of s and s', then M-connecting s to p and p to s' also M-connects s to s' . 

Observe that the three bounding boxes B{o, t), B{o, t'), and B{t, t') have pairwise non-empty in- 
tersections. By the Helly property of axis-parallel d-dimensional boxes, there exists a point min(t, t') 
that simultaneously lies in all three boxes. 

M-connecting min(t, t') with t and t' yields 7r(t, t'), and M-connecting min(t, t') with a yields an 
M-path between any two of the points t,t',o by a repeated application of the above observation. 
This completes the proof. D 
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Proof (of Lemmal^. We double the edges of B and construct a Eulerian cycle C that traverses the 
terminals in T U {o} in some order to,ti, . . . , tn- The length of C is at most 2\\B\\ by construction. 
Now consider the shortcut cycle C in which we connect consecutive terminals tj, tj+i by the M-path 



7r(ti,ti+i) as defined in Lemma 10 we set tn+i := to- Clearly ||C|| < ||C||. We partition C it into 
two halves; Cq = {vr(t2i, t2j+i) \ < i < n/2} and Ci = {vr(t2j+i, ^21+2) | < i < n/2 — 1}. For at 
least one of the two halves, say Co, we have ||Co|| < ||-B||. 

We use Co as a partial solution and recursively M-connect the points in the set T' := {min(f2i, t2i+i^ 



< i < n/2}, which lie in Co (see Lemma 10 ), to the origin by an arboresence A'. Lemma 10 implies 
that the resulting network A = CqL) A' is in fact a feasible RSA solution. The length of A is at 
most IIColl + P'll < \\B\\ + \\A'\\. Note that \T'\ < {\T\ + l)/2. 

To summarize, we have described a procedure that, given the rectilinear cycle C traversing 
terminal set TU {o}, computes a shortcut cycle C, its shorter half Co, and a new point set T' that 
still has to be M-connected to the origin. We refer to this procedure as shortcutting. 

To compute the arboresence A' , observe that C is a rectilinear cycle that traverses the points 
in T' . Shortcutting yield a new cycle C of length at most ||C|| < ||C||,a half Cq no longer than ||-B||, 
which we add to the RSA network, and a new point set T" of cardinality \T"\ < \T'\/2 < (|T| + l)/4, 
which we recursively M-connect to the origin. 

We repeat the shortcutting and recurse. Each iteration halves the number of new points, so the 
process terminates in O(logn) iterations with a single point t. Since min(o,p) = o for any point 



p (see proof of Lemma 10) and our original terminal set T U {0} contained o, we must have that 
t = 0. This shows that the computed solution is feasible. As each iteration adds length at most 
||-B||, we have \\A\\ < [log2ra] • \\B\\. D 

B Running Time Analysis 

We first analyze the running times of the algorithm for d > 2 in Section [3] 

Given an instance R of 0-separated d-dimensional GMMN, the algorithm uses d recursive pro- 
cedures to subdivide the problem into d-separated instances. For j G {0, . . . , d— 1}, let Tj{n) denote 
the running time of the j-th recursive procedure. The j'-th recursive procedure takes a j-separated 
instance R as input and partitions it into two j-separated instances, each of size at most \R\/2, 
and one {j + l)-separated instance of size at most |-R|. The partitioning requires 0{n) steps for 
finding the median of the j-th coordinate value of terminals in R. The two j-separated instances 
are solved recursively and the (j + l)-separated instance is solved with the (j -|- l)-th recursive 
procedure. Let T(i{n) denote the running time to solve a d-separated instance. As pointed out in 
Appendix ffl we can approximate RSA in (i > 2 dimensions by applying (an extension of) the 
algorithm of Rao et al. |RSHS92] to each orthant separately. This requires O(nlogn) time as does 
the original algorithm. Thus we have 

Td{n) = 0(n log n) 

Tj (n) = 2Tj (n/2) + T,+i (n) for j € {0, . . . , d - 1} 

The running time of our overall algorithm is given by ro(n). Solving the recurrences above yields 
ro(n) = 0(n log ^ 11). The improved approximation algorithm of Theorem pi uses Aroras PTAS 
|Aro98] for rectilinear Steiner trees, which worsens the running time substantially but still leads to 
a polynomial running time. 
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Now we analyze the running time of the improved algorithm of Section |4j Stabbing x-separated 
instances can be done with a sweep-line algorithm in 0(n log n) time. The PTAS for RSA requires 
time 0(n"'^'^ logn) for any e with < e < 1. Hence, we have that Ti{n) = O(n^'^logn), and 



To(n) = 2To(re/2) +Ti(n). Solving the recursion yields a running time of To{n) 
the improved algorithm. 



0{n}-l'^ log^ n) for 



C Example Showing the Tightness of Our Analysis 

Observation 1 There are infinitely many instances where the O {log n)- approximation algorithm 
for 2D-GMMN described in Section^has approximation performance Qilogn). 

Proof. We recursively define an arrangement A{n) of n rectangles each of which represents a ter- 
minal pair; the lower left and upper right corner of the rectangle. By a ■ A{n) we denote the 
arrangement A{n) but uniformly scaled in both dimensions so that it fits into an a x a square. Let 
e > be a sufficiently small number. 

The arrangement A{0) is empty. The arrangement A(n) consists of a unit square Sn whose upper 
right vertex is the origin. We add the arrangement bright '■= £■ A{{n— l)/2) and place it in the first 
quadrant at distance e to the origin. Finally, we add the arrangement ^left := (1 — e) • A{{n — l)/2) 
inside the square Sn so that it does not touch the boundary of Sn- See Fig. [5] for an illustration. 
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For each terminal pair 



Fig. 4: Illustration of Lemma 
t,t' , we compute a suitable point min(i, i') and an M- 
path TT{t,t') containing min(f,f'). Adding an arbitrary 
M-path from min(t, t') to o M-connects t and t' to o. 
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Fig. 5: Recursive construction of the arrangement A{n). 
The gray M-path n„ shows an optimum solution. The 
dashed vertical line marks where the algorithm sepa- 
rates. 



Let p{n) denote the cost produced by our algorithm when applied to A{n). Observe that our 
algorithm partitions A{n) into subinstances i?ieft = ^left, Rmid = {Sn}, and i?right = bright- Solving 
the x-separated instance Rmid by our stabbing subroutine costs 1. Let p{n) be the cost of the 
solution to A{n) that our algorithm computes. Recursively solving Rien costs (1 — e) • p{{n — l)/2). 
Recursively solving i?right costs e ■ p{{n — l)/2). Hence, the cost of the solution of our algorithm is 
p{n) > 1 + p{{n — l)/2). This resolves to p{n) = i7(logn). 

Finally, observe that the optimum solution is a single M-path 7r.„ of length 1 -|- 2e going from 
the third to the first quadrant through the origin, see Fig. [5] D 
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